Mean field loops versus quantum anti-crossing nets 
in trapped Bose-Einstein condensates 
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We study a Bose-Einstein condensate trapped in an asymmetric double well potential. Solutions 
of the time-independent Gross-Pitaevskii equation reveal intrinsic loops in the energy (or chemical 
potential) level behavior when the shape of the potential is varied. We investigate the corresponding 
behavior of the quantum (many-body) energy levels. Applying the two-mode approximation to the 
bosonic field operators, we show that the quantum energy levels create an anti-crossing net inside 
the region bounded by the loop of the mean field solution. 
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Experimental realizations of a Bose-Einstein conden- 
sate (BEC) in cooled and trapped atomic gases [1] are 
triggering many efforts to explore new quantum phenom- 
ena on a macroscopic level. The weak interaction be- 
tween the condensate atoms allows for a precise tailoring 
that is not feasible experimentally with other superfluid 
systems such as the ''He-II. In particular, interference 
phenomena have been investigated in the case of a con- 
densate trapped in both double [2] and multiple [3] wells 
potentials, and the possibility of tunneling has resulted 
in an oscillatory (Josephson) atomic flow across adjacent 
traps [4]. A theoretical description of the process car- 
ried out within the mean field (Gross-Pitaevskii equa- 
tion (GPE) [5]) approximation predicts a rich dynamical 
phase diagram [6]. Quantum corrections show that the 
mean field description remains valid for a very long but 
finite period of time up to the moment when the system 
is able to detect tiny differences in energy between the 
multi-particle levels [7-9]. 

The excitation of solitons and vortices in a conden- 
sate by an adiabatic deformation of the trapping poten- 
tial has been recently investigated in [10,11]. However, 
we have found that some levels of the time-independent 
GPE reveal loop-like structures when varying a parame- 
ter of the potential [11]. Solutions of GPE which corre- 
spond to eigenvalues on those loops are typical of non- 
linear equations and do not have a linear (go = in (4)) 
single particle counterpart. This class of GPE solutions 
has been investigated in [12] for a symmetric double well 
potential. The loop-like behavior in the solution of the 
GPE has been also observed for Bloch bands spectrum 
in Ref. [13]. There, the appearance of the loops results in 
the breakdown of the Bloch oscillations due to non-zero 
adiabatic tunneling into the upper band. 

One may wonder about the quantum (multi-particle) 
equivalent of such loops. In this case, indeed, no loops 
are expected because the many-body Hamiltonian is a 
Hermitian operator. It has been proved in [14] that the 
ground state and the associated energy of GPE for repul- 
sive atom-atom interactions (positive scattering length 



tts) are asymptotically exact when the number of parti- 
cles N ^ oo while Nug is kept constant. To the best 
of our knowledge there is no similar theorem concern- 
ing the negative case and the excited eigenstates of 
the GPE. Therefore, it is important to understand how 
those mean-field loops are approached by the eigenen- 
ergies of the quantum Hamiltonian which, we remark, 
cannot form any loop. 

The aim of the present paper is to investigate quan- 
tum multiparticle energy levels of an asymmetric double 
well potential in the regime where the GPE reveals loop- 
like structures. The structures can be observed even for 
a very weak interaction between particles, so we use a 
two-mode approximation (TMA) [7,6,8,9] in the quan- 
tum calculations. The results so obtained are compared 
with the GPE predictions, calculated both in the TMA 
and with a numerical integration in a ID geometry. 

We consider a condensate trapped in a harmonic cigar- 
like potential with an additional Gaussian perturbation. 
Freezing the transversal degrees of freedom tightly con- 
fining the condensate in the same direction, the potential 
(in units of a harmonic oscillator corresponding to the 
longitudinal direction) reads 
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Within a certain window of parameters, the potential (1) 
has a double well shape [10,11]. The many-body Hamil- 
tonian governing the full A'^-particle problem [15] is. 



H 



dx 



(2) 



where Hq = —5^^ + V{x). H can be approximated ex- 
panding the field operators in terms of the local modes, 
ip ~ Ui{x)ci + U2{x)c2. Uj{x) are localized solutions 
of the Schrodinger equation solved for each well of the 
potential separately (see below). This, provided that 
J ul{x)u2{x)dx « 0, yields 
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(3) 



where Ej = J u*{x)HoUj{x)dx, = 2/ ul{x)HoU2{x)dx 
and the effective mode volume of each well = 

IMx^dx. 

The Bose operators cj , that annihilate a particle at the 
respective well of the potential, obey the usual commuta- 
tion relation [cj, c^] = Sjk- To define the modes Uj{x), we 
have separated each well from the potential (1) as shown 
in Fig. 1 [17]. Given the modes one can immediately 
calculate the parameters of the Hamiltonian (3). 
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FIG. 1. (a): plot of tfie potential (1) for Uo = 6.4, a = 0.5 
and xo = —3.7. (b): single well potential obtained by tak- 
ing the left part of the potential plotted in panel (a) up to 
the point corresponding to the top of the barrier and in the 
remaining space the potential is kept constant. The lowest 
eigenstate of such a single well potential constitutes the mode 
ui{x). In an analogous way the mode U2{x) of the right po- 
tential well is defined. 



The solutions of the time-independent GPE, 

H(np{x) + NgQ\Lp{x)\^ip{x) = iJ.(p{x), 



(4) 



where |¥'(a;)p is normalized to unity and N denotes the 
number of particles in a BEG, give the chemical poten- 
tial levels /U of the condensate [5]. For a very small 
nonlinearity parameter, Ngo, changing the parameter of 
the potential xo (for Uq = 6.4 and a = 0.5), a nar- 
row avoided crossing between the ground and first ex- 
cited levels occurs around .tq = —3.6 [10]. However, 
when Ngo increases the usual avoided crossing turns into 
a more complicated structure where the excited level 
reveals a loop-like character, sec Fig. 2. To compare 
this single-particle description of the system with the 
quantum multi-particle picture we have diagonalized the 
Hamiltonian (3) for = 10 and = 50 adjusting go so 
that Ngo = 1 in both cases. Figure 3 shows that quan- 
tum energy levels undergo a net of anti-crossings and the 
GPE solutions provide an envelope of this anti-crossings 
structure. The lowest quantum energy level corresponds 
to the ground state of the BEG while the highest one 



corresponds to all particles occupying the upper level of 
the trap. Both states can be calculated by solving the 
GPE equation. However, in the middle of the figure not 
only the lowest and highest energy states can be repre- 
sented by a mean field wavefunction which corresponds 
to all atoms in the same one-particle state. Indeed, the 
existence of the loop indicates that there are two more of 
such states. 
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FIG. 2. Energy levels of a condensate per particle as a 
function of xo, for Ngo = 1. The other potential param- 
eters are: Uo = 6.4 and cr = 0.5. Thick lines, in both 
panels, correspond to the solutions of the GPE (4), i.e. 
E/N = n - {Ngo/2) J \ip{x)fdx, while lines are the results 
of the quantum two-mode approximation (3) for number of 
particles N = 10 (a) and iV = 50 (b). 



As expected, also in the attractive case we observe a 
loop-like structure (see Fig. 4) which, however, is on the 
ground state energy level. In the two-mode approxima- 
tion, the conditions for the formation of loops can be 
easily calculated analytically (see below). 




FIG. 3. Magnified anti-crossings structure of Fig. 2b. 
Anti-crossings are so tiny that they become visible in the in- 
sert (in the circle). 
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Switching to the quantum two-mode approximation we 
found that the energy levels reveal (similarly as for the 
repulsive interaction) an anti-crossing structure located 
within the loop of the mean field results (Fig. 4). 
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FIG. 4. The same as in Fig. 2 but for Ngo = -1. 

Let us now carry out the two-mode approximation also 
in the mean field description. Replacing the ansatz wave- 
function ip{x) = {ciui{x) + C2U2{x)) / \/N in (4) we get a 
pair of coupled equations, 



C2^2+Cl§ + |C2|2c2f =C2M ' 



(5) 



which (with the additional constraint |cip -|- |c2p = A'') 
allow us to calculate the coefficients Cj and the chemical 
potential. In the range of the parameter xq where there 
is no loop of the GPE, we get only two solutions of (5), 
while in the inner region there are four of them form- 
ing a loop. In Fig. 5 we compare the classical two-mode 
approximation with the results of the GPE both for re- 
pulsive and attractive interaction between particles. The 
results of the classical two-mode approximation are in a 
good agreement with the full GPE solutions and, there- 
fore, we expect that the quantum two-mode will also be 
a good approximation of the full quantum multi-particle 
problem. We note that the expression for the energy of 
the A^-particle condensate in the classical two-mode ap- 
proximation, 



H = Ei\cif + E2\C2f + ^{ci4 + clC2) 
, 90 , |4 , 90 , |4 



(6) 



is of the same form as the corresponding quantum Hamil- 
tonian (3), but with cj being c- numbers. 

By increasing the number of atoms in the BEG keeping 
Nqq fixed, the values of the extremal quantum energy 
levels reproduce the classical two-mode approximation 
results more and more accurately. 
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FIG. 5. Comparison of the two-mode approximation ap- 
plied to the GPE (solid lines) with the solutions of the full 
GPE (dashed lines). Panel (a) corresponds to Ngo = 1, panel 
(b) to Ngo = -1. 



In [10] we have suggested that starting with a conden- 
sate in the ground state of a harmonic trap and sweeping 
the trap with a laser beam (which actually leads to the 
potential (1) where xq changes from some negative value 
to zero) can be used to excite a BEG solitary wave by 
means of a diabatic passage through an avoided crossing. 
In the mean field description, the presence of a loop does 
not break the process because; the ground state wave- 
function on one side of such a strange avoided crossing 
structure (i.e. a loop) is basically the same as the wave- 
function of the first excited state on the other side so the 
diabatic scenario applies [11]. Switching to the quan- 
tum description we see that one avoided crossing of GPE 
levels, which has to be passed by each one of the N par- 
ticles, is just replaced by many smaller ones in the multi- 
particle level structure. We expect that when the laser 
beam is swept with appropriate velocity it will transfer 
the atoms from the ground energy level on the left hand 
side of Fig. 2 (or Fig. 4) to the highest energy levels on 
the right hand side of Fig. 2 (or Fig. 4). To test such a 
hypothesis we have carried out quantum time dependent 
calculations within the TMA starting with the system in 
a ground state of the pure harmonic trap [18]. In Fig. 6 
we show the results of a laser sweep, i.e. a change of xq 
from —5 to according to xo{t) = —5 + 0.05t. In the 
case of the noninteracting particle system (go = 0), if the 
excitation probability of a particle from the ground state 
to the first excited state in the single particle description 
is less than 1 (although close to unity), the correspond- 
ing multiparticle picture reveals population of a bunch of 
the eigenstates as shown in Fig. 6a [10] . Switching on the 
interaction the effect of the laser sweep is unchanged — 
similarly as in the noninteracting particle case the high- 
est energy eigenstates of the system become populated, 
see Fig. 6b and Fig. 6c. Of course, this is a simplified 
picture in the two-mode approximation that is valid only 
for a weak interaction between atoms. For stronger in- 
teractions further investigations are necessary. 
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FIG. 6. Overlaps |(n|?/>)p of the wavcfunction on the en- 
ergy eigenstatcs of the A'' = 50 particle system in the har- 
monic trap at the end of the laser sweep. Starting with the 
ground state of the harmonic trap, i.e. |V'(0)) = |n = 0), we 
change the potential (1) according to xo{t) = a;o(0) + 0.05t 
with a;o(0) — —5. At the end of the sweep the highest en- 
ergy eigenstates become populated — panel (a) corresponds 
to the noninteracting particle case (i.e. go =0), panel (b) to 
goN = 1 and (c) to goN = -1. 

In Ref. [13] a loop-like behavior was observed in solu- 
tions of the GPE for a BEG in a periodic lattice, implying 
breakdown of Bloch oscillations and adiabatic tunneling 
into the upper band. The approximation method per- 
formed here can be easily applied in that case. This al- 
lows one to identify which multi-particle eigenstates are 
populated in the tunneling process. 

To summarize, we have considered the behavior of 
the energy levels of an asymmetric double well potential 
when the shape of the potential is changed. Using the 
two-mode approximation in the quantum multi-particle 
Hamiltonian we have found that the energy levels un- 
dergo a net of avoided crossings. Solutions of the Gross- 
Pitacvskii equation reveal a loop structure which pro- 
vides an envelope to the quantum avoided-crossings net. 
We finally note that since the two-mode approximation, 
applied to the mean field description, is in a very good 
agreement with the full numerical solution of the Gross- 
Pitaevskii equation, the quantum two-mode approxima- 
tion can also be considered a good approximation of the 
full quantum multi-particle problem at least in the regime 
of weak interaction between atoms stiidicd here. 

Discussions with Bogdan Damski, Jacek Dziarmaga 
and Kuba Zakrzewski are acknowledged. KS is grate- 
ful to Wojcicch H. Zurck for the hospitality during the 
preparation of this work. Support of KBN under project 
5 P03B 088 21 (KS) and of by the U.S. DOE are acknowl- 
edged. 



[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. 

Wicman, and E. A. CorneU, Science 269, 198 (1995); 

K. B. Davis, M.-O. Mewcs, M. R. Andrews, N. J. 

van Druten, D. S. Durfee, D. M. Kurn, and W. Ket- 

terle, Phys. Rev. Lett. 75, 3969 (1995); C. C. Bradley, 

C. A. Sackett, J. J. Tollett, and R. G. Hulet, Phys. Rev. 

Lett. 75, 1687 (1995); 
[2] M. R. Andrews, C. G. Townsend, H.-J. Miesner, D. S. 

Durfee, D. M. Kurn, and W. Ketterle, Science 275, 637 

(1997) . 

[3] B. P. Anderson and M. A. Kasevich, Science 282, 1686 

(1998) ; C. Orzel, A. K. Tuchman, M. L. Fenselau, M. 
Yasuda, and M. A. Kasevich, Science 291, 2386 (2001). 

[4] F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. 

Minardi, A. Trombettoni, A. Smerzi, and M. Inguscio, 

Science 293, 843 (2001) 
[5] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 

Rev. Mod. Phys. 71, 463 (1999). 
[6] A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy, 

Phys. Rev. Lett. 79, 4950 (1997); S. Raghavan, A. 

Smerzi, S. Fantoni, and S. R. Shenoy, Phys. Rev. A 59, 

620 (1999) 

[7] G. J. Milburn, J. Corney, E. M. Wright, D. F. Walls, 

Phys. Rev. A 55, 4318 (1997). 
[8] S. Raghavan, A. Smerzi, and V. M. Kcnkre, Phys. Rev. 

A 60, R1787 (1999). 
[9] M. Holthaus, Phys. Rev. A 64, 011601(R) (2001). 
[10] Z. P. Karkuszewski, K. Sacha, and J. Zakrzewski, Phys. 

Rev. A 63, 061601(R) (2001). 
[11] B. Damski, Z. P. Karkuszewski, K. Sacha, and J. Za- 
krzewski, Phys. Rev. A 65, 013604 (2002). 
[12] R. D'Agosta and C. Presilla, Phys. Rev. A 65, 043609 
(2002). 

[13] Biao Wu and Qian Niu, Phys. Rev. A 61, 023402 (2000). 

[14] E. H. Lieb, R. Seiringer, and J. Yngvason, Phys. Rev. A 
61, 043602 (2000). 

[15] In the harmonic oscillator units, go = g/{hu!\\a\\) where 
the coefficient g = ATrh^as/imna^) is renormalized by 
the transverse condensate length a± in reducing the di- 
mensionality of the problem from 3D to ID [16]. 

[16] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

[17] The modes may be defined in many ways, e.g. as states of 
potentials obtained by expanding Eq. (1) around a cen- 
ter of each well up to the second order (which, however, 
is a poor approximation for the Gaussian-like well which 
we deal with in our problem) or as symmetric and anti- 
symmetric eigenstates combinations. Such states can be 
calculated as solutions of the Gross-Pitaevskii equation 
(providing, in general, a more accurate two-mode approx- 
imation [6]). In the present paper, however, the localized 
wave-functions have been calculated by neglecting the (in 
our case very small) non-linearity. 

[18] In the time dependent calculations we have defined 
modes as the two lowest eigenstates of the potential (1). 
It allows us to describe the system not only when the 
potential (1) reveals double well structure but in the full 
range of the potential parameter xo- 



4 



